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ABSTRACT 

A mathematical model to represent hysteresis in dynamic 
systems was proposed in a recent paper by Chua and Stromsmoe [8]. In 
the case of iron-core inductors and transformers the model has an 
extremely simple circuit interpretation, namely a parallel combination 
of a nonlinear resistor and nonlinear inductor and ideal trans- 
formers. In the previous work several experimental results were 
compared to those predicted by the model. All examples were limited 
to circuits containing either an iron-core inductor or an iron- 
core transformer. 

In the present work the same model is used to examine 
circuits containing two iron-core transformers. The specific 
circuits considered are the parallel and series connected non- 
feedback magnetic amplifiers. In both circuits the magnetizing 
current is appreciable in comparison to the total current. Thus the 
model must accurately account for both large nonlinear effects and 
appreciable hysteretic losses. 

The magnetic amplifier is operated under a wide range 
of conditions, the simulated and measured results are compared. The 
close agreement between these indicate that the model is remarkably 


good in predicting the performance of the magnetic amplifier. 
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CHAPTER 1 


Tet Tntroduc tion 


Hysteresis is a phenomenon that is observed in many 
materials and processes. Hysteresis is observed in the flux versus 
magnetization force relationship of ferromagnetic and ferrimagnetic 
materials [1-3], in dielectrics [4], and in the stress-strain re- 
lationship of materials that undergo plastic deformation. The 
human respiratory system is another example of a process where 
hysteresis occurs [5]. In order to analyze systems, in particular 
circuits, containing hysteretic elements, it is necessary to have 
a mathematical model which can replace the hysteretic elements and 
closely simulates the hysteresis phenomenon. 

To solve this modeling problem several mathematical 
models have been reported. These mainly follow two basic approaches, 
the first approach [6] is a physical one based on the physical 
principles governing the behavior of the hysteretic elements, and 
the second approach [7] is based on postulating a suitable mathematical 
representation which can be shown to exhibit hysteresis phenomenon. 

The first approach results in mathematical models in the form of 
partial differential equations, and even with the help of a computer 
the analysis of nonlinear systems or eee containing the hysteretic 
elements described by these differential equations will be difficult. 


The second approach results in nonlinear mathematical models. 
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The model used in the present work is based on this approach [8]. 
The establishment of the accuracy and validity of these models in 
simulating hysteretic elements is not an easy process. For 

linear models there is a valid representation theorem [9] stating 
that if a linear model is found to give an exact response, such 

as impulse response, corresponding to a given test signal, then the 
model will give the correct response to all other test signals. 
This representation theorem is not valid for nonlinear models since 
its proof makes use of the principle of superposition. Without a 
representation theorem, analogous to that of linear systems, it is 
difficult to establish the validity of nonlinear models. This is 
the main difficulty inherent in the second approach where the 
establishment of the validity of the nonlinear model requires an 
infinite set of measurements corresponding to all possible excitation 
signals. 

Actually the validity of any of the models based on 
this second approach is first established qualitatively by showing 
that the postulated model exhibits the same significant properties 
and features of the hysteretic elements, such as expansion of the 
loop with increase in frequency and the presence of minor loops. 
The second step is to construct the functions and parameters of the 
mathematical representation so that the resulting model would yield 
realistic response to one or more test signals. The model is then 


assumed to yield reasonably accurate responses to other excitations. 
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The exactness of the model in representing the hysteretic elements 
depends. on its accuracy in predicting the correct solutions when 
imbedded in a large number of systems under a variety of excitations. 
Therefore the ultimate verdict on the usefulness of the model will 
not be known until sufficient data is gathered concerning its 
ability to simulate the correct solutions to a large number of 


practical problems. 
1.2 Objective of this study 


A mathematical representation based on the second approach 
was reported recently [8]. This model is a simple one and is very 
useful for analysis of practical circuits containing hysteretic 
elements. A major advantage of this mathematical representation 
is that all its functions can be determined from a given measurement. 
This mathematical model is an important and useful one especially in 
the case of iron-core inductors. For such inductors the mathematical 
representation is found to be equivalent to a simple lumped circuit 
model. This lumped circuit model makes it possible to analyze 
circuits containing iron-core inductors or iron-core transformers by 
simple computer programs. 

This mathematical model was proved qualitatively to exhibit 
most of the significant properties and features of the hysteretic 
elements. As a step in proving how accurate this model is, it was 
shown [8] that this model yields realistic results when used to represent 


the hysteretic elements in a push-pull tunnel diode relaxation 
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oscillator which use a saturable reactor and tunnel diodes. The 
model was found to give a very good prediction concerning the frequency, 
the waveform and amplitude of the generated signal. 

As a second step in eheeiine the exactness of this 
mathematical model, the present work is intended to check its ability 
and usefulness in analyzing another practical circuit which is the 
Magnetic amplifier. In the case of the push-pull tunnel diode re- 
laxation oscillator there was only one iron-core transformer while 
in the case of the magnetic amplifier there are two iron-core 
transformers. The magnetic amplifier is analyzed using the pre- 
described model to represent the two iron-core transformers. The 
waveforms and amplitudes of the different voltages in this practical 
circuit are computed numerically. These are then compared to the 
actual experimental waveforms and amplitudes of the corresponding 


voltages, the agreement is found to be excellent. 
1.3 Magnetic amplifiers 


Magnetic amplifiers are devices which use a small 
applied input signal to control the voltage or power output. The 
input signal is at a power level much lower than that of the output 
power. The magnetic amplifier consists of a combination of saturable 
reactors, rectifiers, resistors and transformers. The operation of 
the magnetic amplifier is based on the characteristics of the single 
saturable reactor. A saturable reactor has a ferromagnetic core and 


two windings, one of these is the d.c. or control winding N. and the 
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other is the a.c or gate winding Sd For optimum performance of 
the magnetic amplifier, the core material of its saturable reactors 


should have the following magnetic properties:- 


1. Hysteresis and eddy-current losses a minimum. 

2. High saturation flux density, to obtain large power 
handling capacity of a given weight of core material. 

3. The general shape of the hysteresis loop as close as 
possible to a thin rectangular one. 

4. Stability of magnetic characteristics under changing 


temperature, mechanical strain and shock conditions. 


The last property ensures that the characteristics shall 
not change under varying operating conditions. The first three 
properties ensure that a small change in the control current (input), 
flowing through the control winding, varies the magnetization state 
of the ferromagnetic core over wide limits. This results in a wide 
range of variation of the inductive impedance of the gate winding. 
Thus the impedance of the gate winding is dependent on and controlled 
by the value of the d.c. control current. This principle was used 
to construct different types of magnetic amplifiers. In all these 
types the magnitude of the control signal controls the amplitude and 
the waveform of the voltage across, and hence the power delivered 
to,the load. This controllability together with the high ratio of the 


output power to the load to the input power to the control circuit 
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provide the desired controllable high gain in magnetic amplifiers. 

Like other types of amplifiers the gain of the magnetic 
amplifier increases when positive internal (intrinsic) or external 
(extrinsic) feedback is introduced. Generally the single-stage 
and multi-stage magnetic amplifiers used in commercial and industrial 
applications make use of this positive feedback to allow for higher 
power amplification [10]. However, it must be noted that the 
operation of these feedback magnetic amplifiers is based on the 
characteristics of the non-feedback type. In addition, it must be 
noted that the presence of the many applications of the feedback 
pacnetic amplifiers does not imply that those without feedback do 
not serve a useful purpose. Despite the lack of high power gain 
of the latter, they are sefatene simple and have great versatility 
and do have certain important and useful applications [10-11]. 
Actually a common factor amongst all types of magnetic amplifiers 
is that they have extremely high reliability, long life, ruggedness, 
no warm-up time, high efficiency, a minimum maintainance problem 
and are capable of high temperature operation. 

The present work deals only with the parallel and series 
connected non-feedback type of magnetic amplifiers. As their names 
imply, the series-connected type has its gate windings connected in 
series while the parallel-connected type has its gate windings 
connected in parallel. The parallel-connected magnetic amplifier 


has a much slower response to control signals than that of the series- 
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connected magnetic amplifier [10]. The parallel-connected magnetic 
amplifier is therefore used in high-power applications where the 
speed of the amplifier response to changes in control signal is 
not of prime concern, while the series-connected magnetic amplifier 
is usually used in applications where nominal power is required and 


where maximum control speed is required. 
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CHAPTER 2 
THE MATHEMATICAL MODEL 
2.1 Mathematical representation 


Hysteresis, as stated before, is a phenomenon exhibited 
by many devices. The present work deals mainly with the hysteresis 
observed in magnetic materials, for example iron-core inductors or 
transformers. For magnetic materials the two variables exhibiting 
the hysteresis phenomenon are the flux linkage A(t) and the wind- 
ing current i(t) causing this flux linkage to exist. The mathematical - 
model, whose accuracy in representing magnetic materials is to be 
examined in this work, is a dynamic model. This model relates the 


two variables \(t) and i(t) through the differential equation, 


Swedes eQuli ftiyutet (ad) (2.1) 


This differential equation completely defines this mathematical 
model. The two functions f(-) and g(+) are real - valued functions 
defined on the real line R and the symbol "o" denotes the composition 
operation. The detailed properties of this mathematical model, 
together with the proofs are in reference [8]. The proving of these 
properties depends on the fact that the two functions f(-) and g(-) 


should satisfy the following conditions, 
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(i) They are strictly monotonically increasing functions. 
(ii) Both are differentiable onto functions. 
(i.43.) The two functions have non-zero slopes throughout the 


entire real line. 


(iv) They should satisfy the property g(0) = £(0) = 0. 


The function g(-) is called the dissipation function since it can be 
shown to be responsible for the energy dissipation in the element 

while f(*) is called the restoring function since it is responsible 
for the energy storage in the element. This mathematical model can 
represent any of the hysteretic magnetic materials. In the present 


work only iron-core transformers are to be replaced by such a model. 
2.2 Iron-core transformer 


The mathematical model defined by equation (2.1) has 
Many advantages, the most important one is that it yields a lumped 
circuit model that simulates the hysteretic behavior of an iron- 
core inductor. This lumped circuit model is shown in Fig. 2.1. 
Unlike the linear inductors and resistors, it is necessary for 
nonlinear resistors and inductors to indicate the positive direction 
of the current flowing through each. Thus each of the two non- 
linear elements in Fig. 2.1 is represented by a rectangular symbol 
having one of its ends darkened. This dark end is the negative end 
of the nonlinear element and the positive current flows from the 


positive end to the negative end. Inside the rectangular symbols 
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the conventional representations of inductors and resistors are used. 
The lumped circuit shown in Fig. 2.1 consists of a nonlinear inductor 
"L" in parallel with a nonlinear resistor "R". The nonlinear 
resistor "R'" is defined by the relation between its current in and 
the voltage across its terminals VR: The characteristic equation of 


this resistor is, 
it = = = = 7 
v(t) = v. Vp = 8 (42) (252) 


The nonlinear inductor "L" is defined by the relation between its 
current i and the flux id linking the iron-core inductor. The 


characteristic equation of this inductor is, 


i, = £Q) (2.3) 

It is clear from this lumped circuit model that the function g(°) 

is responsible for the energy dissipation in the iron-core inductor, 
while the function f(*) is responsible for the energy storage. This 
circuit model of the iron-core inductor is used to obtain the circuit 
model for an iron-core two winding transformer. Thus the trans- 

former will be completely represented by the circuit model shown in 
Fig. 2.2. The nonlinear resistor and inductor are defined by the 
characteristic equations (2.2) and (2.3) respectively, where v(t) 


is the voltage across one of the two windings of the transformer and 


» 


chew oe ex0snieod ‘pale 2 sugar i 
sodeuibat sponblaco 5 Qo aaatenos. hey 
ent Ln oo Hey Hegde: 


bare b Inbrsu3 eit Baas ried ote oe we - ee 


a 
to ‘aol isups sivaiaayen seo ost 


‘tes ‘sodoubat hay 20 sot teups 
z pie oe a 


(€.8) 


cs a. Olaga0 ‘ota Soa fobon teszae Saget ere 
.rosoubn e1d-noxt at. ab nobisgiee 
shit } -auering at oa: 3 


stupaken eda, stele : bi 


=A 


ARM 


PIG... 261 


Lumped circuit model for an iron-core inductor 


FIG, 2.2 


Lumped circuit model for an iron-core transformer 
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A(t) is the flux linking this winding. This is the circuit model 


that is used in the analysis of the magnetic amplifiers. 
2.3 Determination of f(:) and g(-) 


Since the model was found to exhibit most of the important 
properties observed in practice [8]. Hence if f(-) and g(-) are 
identified to simulate the hysteresis loop of the iron-core inductor 
under consideration, then it is reasonable to assume that the model 
will give realistic results under any type of excitation. The type 
of excitation is dependent on the system in which the iron-core 
inductor or transformer is imbedded. To identify f(-) and g(-) a 
pair of waveforms i(t), A(t) must be measured. The nonlinearity 
of the model implies that there is no apparent advantage of one 
set of measurement over the other. Thus, the obvious choice is to 
use a set of measurements that can easily lead to the identification 
of the nonlinear functions f(-) and g(-). 

Since the iron cores used in magnetic amplifier have 
symmetric hysteresis loops, the two functions f(-) and g(+) are odd 
functions [8]. Assuming a cosine like variation of A(t), a sinusoidal 
test signal being used to produce the hysteresis loop, then A(t) is 
an even function of time while its derivative \'(t) = < ist an*odd 


function of time. Since f(-) and g(-) are odd functions, 


EO) ac) &) (2.4) 
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Thus f(\(t) ) is an even function of time while aa O:'.Cb))) sts.an 


odd function of time. If equation (2.1) is rearranged in the form 


i(t) =g 7 GMt)) + £ C(t) (2.5) 


and the current i(t) is expressed as the sum of an even component 


i (t) and an odd component i (t) where, 


ote i(t) + i(-t) 
e 


Z 
(2.6) 
, be V(t), =e (ee) 
i (t) 7 2 
Then, equations (2.4), (2.5) and (2.6) show that 

i(t) = f(r (t)) (2-7) 
° 1 1 
1 kt Ome BE RINE CC ):) (2.8) 


For X(t) a cosine function of time, there are two instants ty and t, 


in each cycle at which 


A = A(t) = (ty) (2.9) 
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ee AN(t,) = A" (ty) CZ 
Equations (2.7) to (2.10) give 

i(t,) = £(A(t,)) = £(A(t,)) = i (t,) el. C2. 

i (t,) =e Q'(t,)) = -g @'(t,)) = -i.(t,) (2 

Ope 2 1 Ouse: ; 
But 

i(t,) = i (t,) + i(t,) (2% 

i(t,) = ii (t,) + i(t,) (2. 
From (2.11) and (2.12) equation (2.14) becomes 

i(t,) = - i (ty) + i(t,) (2. 


Thus trom. €2.13)--and, (Zeno) (A) and ee) can be related to i(t,) 


and i(t,) as follows, 


£0) = FA(t,)) = a.ce,) = i (ty) 4 i (2) =i. (2 
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positive or negative sign is dependent on the value of A'. Since 


In equation (2.17) d°= 
g(-)is a strictly monotonically increasing function and g(0) = 0 then 
eee oa eae? ee SiCincity 
g (A') = -d for} A' < 0 


Equations (2.16) and (2.17) show that through a set of measurements 
and using these two equations the two functions f(-) and g(*) can be 
completely defined. Thus to determine f(+) and g(:) the hysteresis 
loop is obtained using the sinusoidal test signal. For each value 

of A} a horizontal line is drawn intersecting the hysteresis loop at 
two points corresponding to currents i(t,), i(t,), ty and ty being the 
two instants at which A(t) = (to) and equal to the specific value 

of X} considered. The value of the current i, is calculated as the 
average of i(t,) and i(t,) while d is half the magnitude of the 
difference between them. If this process is repeated, for values 

of X spanned by the hysteresis loop, the dashed line in Fig. 2.3 is 
obtained which lies midway between the locus of i(t,) and i(t,). This 
dashed line represents the relation meee i), and the nonlinear 


function f(-) is defined. For each value of X considered the value 
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BIG. 72.3 
Construction procedure for determining the restoring function f(.-)and 


the dissipation function g(-) 


i] 


of A' is known and the corresponding value of d is already obtained. 
Thus the relation between d and ' can be plotted which completely 
defines the nonlinear function g(-) as in equation (2.17). Once 
f(-) and g(-) are defined, this completely defines the character- 
istics of the nonlinear inductor and resistor in the model of the 
iron-core transformer. 


2e4 Model for the magnetic amplifier cores 


The parallel and series connected magnetic amplifiers 
analyzed in the present work both use two saturable reactors having a common 
control winding. The circuit model for these two saturable reactors 
with their three windings is shown in Fig. 2.4, where the transformers 
shown are ideal transformers, 
Since the two saturable reactors are identical the 
two nonlinear inductors L, and L, are the same and the two nonlinear 


‘} 2 


resistors Ry and Ry are the same. To determine the characteristics 
of these nonlinear resistors and inductors the dissipation function 
g(-) and the restoring function f(*) should be identified. The 
hysteresis loop was experimentally obtained. Then following the 
steps described in the previous section f(+) and g(-) were 
obtained. The steps predescribed were used to calculate data 
points on the nonlinear characteristics of f(-) and g(-). In the 


present work these steps were carried out at eight values of i. 


These eight values of X were chosen so that the whole range of 
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variation of X is covered, from negative saturation up to positive 
Saturation. Thus eight data points were obtained on the nonlinear 
characteristics of the nonlinear inductor and the nonlinear resistor, 
Table 2.1: In order to calculate intermediate points some inter- 
polation method must be used [12]. Due to certain reasons discussed 
later in Chapter 3, it was found that the piecewise linearization 
technique has to be used to obtain intermediate points, between the 
eight data points, on the characteristic curves of the nonlinear in- 
ductor and resistor. Thus the nonlinear characteristic curves were 
replaced by piecewise linear characteristic. These piecewise linear 
characteristics were then used to calculate and plot the hysteresis 
loop. The resulting hysteresis loop is shown in Fig. 2.5 which 
clearly agrees with the actual experimental hysteresis loop shown 

in Fig. 2.6. In both figures the vertical scale is 0.01 Weber/division 


and the horizontal scale is 0.03 ampere / division. 
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FIG. 2.5 


Hysteresis loop computed using the lumped circuit model 


FIG. 2.6 


Experimental hysteresis loop 
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CHAPTER 3 
PARALLEL-CONNECTED MAGNETIC AMPLIFIER 


3.1 Theory of operation 


The circuit diagram of a parallel-connected non-feedback 
magnetic amplifier is shown in Fig. 3.1. This is an elementary but 
practical circuit [10-11-13]. As mentioned previously in Chapter 2 
the two saturable reactors have a single control winding of Nc turns. 
This control winding is common to the gate windings of the two 
a aN reactors. The two gate windings, each having turns, are 
connected in parallel in such a way that a voltage applied to 
their parallel connection will, theoretically if Re Re, ; 
induce no a.c. current in the control circuit. The two, parallel- 
connected, gate windings are in series with a resistive load Ry 
and in series with an a.c. voltage source es: The magnetic amp- 
lifier output is the a.c. voltage across the load Ro: liane. 
output is required, a bridge rectifier can be used. The rectifiers 
ane affect the operation or the analysis of the magnetic amplifier 
and are merely considered as a part of the doad. The control winding 
No is in series with the d.c. control voltage EQ + and in series with 


'a large resistor. The resistance Ro is the sum of the d.c. control 


winding resistance plus that of the series large resistor. The large 


value of Ra helps to reduce the value of the a.c. current induced in 
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the control circuit by transformer action. 

For Eq=0 the control winding contributes m d.c. flux 
to the cores and the total flux in each of the cores is produced by 
the a.c. current only and is much less than the value of the flux 
at saturation. Therefore, the permeabilities of the cores and hence 
the input impedances of the two gate windings are high compared to 
R. and only a very small portion of the a.c. voltage eg appears 
across R For Eo + O, there is some d.c. flux linkage in each 


of the two windings due to the d.c. control current i There is 


Cc: 
also some a.c. flux linkage in each of the two windings due to the 
a.c. current in the load circuit. In each half of the cycle one 

of the two gate windings will have the d.c. flux linkage in it 

due to in add to that due to the a.c. load current and the total 
flux linkage in this winding is driven towards the saturation value. 
This winding will have a lower inductive impedance than when Eo was 
zero. In the other winding the two fluxes will oppose and its 

total flux will still be away from saturation value and this winding 
will still have a high inductive impedance. Actually the low 
impedance gate takes most of the load current resulting in further 
reduction of its impedance. Since the two windings are in parallel, 
the input impedance of their parallel connection during this first 
half cycle will be lower than their parallel impedance when Eo was 


zero. In the next half cycle each of the two windings replaces the 


other in action. Thus the input impedance of the two parallel-connected 


| a hee Rate OERo8 tearsotedss co 13 uD" 2 
stink 2 i is peandtxan0> gotbndw toxdmoo ads nye a 


we Besbborg wt @s%0> ea to tose ar ie eres, ie ~ | 
ott ot Yo gutsy srt, ad eee! fous at bag: ry 1a 


ea a if ‘spas | baé BG2209 Sd to sobsiliddomssg sa otgigaaitt 


69 bataq@os dtd SiS i ia ote ows oda 25 5 


| atseqqs 8 Sasa tov Vo. 5 “eis 0. AOEII0q Lisme aay & a, 


dows mi egedail xl so. a Soe ‘ak ategy ° i a aot 


: at sted /.gh Inez. Loxziaos 198 asfd od oub egctbinw 
ah phot odd o2 ub egaibale OWN a3 to fists TE sgerot xult ean 
Vee 7 : S80 efsy> ods ic tad dn 289 ni stupas) aah sity. arr | 
a iA ji at gaiks suLt 1960 ‘aif avd ite oaitt bate e009 hs 
Pa | "faded sits bas soot pact | ‘Ske ad jos “ay dass, 93 bbe “0a 
Dae vee agkianae edt ebasvnd aavirb et igure katy) ee oheanel 
is . ‘, aw ne nadia eed Men ag evboubey hep ‘e as 
; a 


wot oot wilevaon aii ening datit 6 hen 


tee ak gntsiugex sneazuD bebe eile ae Jeon’ gedsi aaa | 

= brag ai 930 eantbaty ows re ia ai cas al 
tae tits ate robssen | 

La . . Nite cists ae aban. ad Di ¢ 

coat ayatbatw ows oda “20. fas oie Ase oxen Por 

ee 


| satiate ows ota — ) sonebaga! juga sia eyitT nents 
= .*) 


- Nol = ae ¢ 
& 7) , Co 


25 


gate windings remains very low over the whole cycle at a value 
dependent on the magnitude of Eo: This causes an increase of the 
a.c. load voltage across Ro and hence an increase of the a.c. power 
delivered to it to values dependent on Eas 
3.2 Analysis 

The first step in the analysis of this circuit is to 
replace the two iron-core transformers shown in Fig. 3.1 with their 
equivalent nonlinear lumped circuit model shown in Fig. 2.4. The 
characteristics of the nonlinear elements of this lumped circuit 
were obtained in the previous chapter. Thus the parallel-connected 
Magnetic amplifier is represented as in Fig. 3.2. All the trans- 
formers shown in Fig. 3.2 are ideal transformers. The equivalent 
circuit of the parallel-connected magnetic amplifier is simplified 
by refering all the circuit elements to the load side with the 
necessary variation in the parameters, Fig. 3.3. Rearranging the 
equivalent circuit and converting all voltage sources to current 


sources, the equivalent circuit takes the form shown in Fig. 3.4 


where, 
vi = Vv. (37 0) 


vy is the voltage across the nonlinear inductor Ly due to variation 
dr 


x 
in the flux A, linking the first gate winding, y, = aes 3 


i 


V. = -V (52) 
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V5 is the voltage across the nonlinear inductor L, due to variation 
di. 


in the flux ro linking the second gate winding, Vas are 


ae i pe 
E' = ee ° ee e 
Eo No 90 ° Since ¥ 79 12 this case. 
C 
Z 
Noy” OR 

R' =R He | JS 

C No 400 


e 


4 Es Sin wt, the frequency of this a.c. source is 60 Hz. 


The currents i. and i are those through the nonlinear inductors 
2 YM 


Ly and L, respectively while i 


9 and i 


are the currents through 


R R 


1 
the nonlinear resistors R, and R, respectively. 


With reference to Fig. 3.4, applying Kirchhoff current 


law to nodes c, d and a gives 


| 
ey it at ny ik 
— -i aid =V a ey ar Gis> 
Ro R, Ly 6 Ra e a oh, d Ro 
En i 1 ees: 
C 
E er at ares Famili Ee 4 
Reece See . (i | Ve Eee d\n bac aR | (3.4) 
C 2 2 C g Cc g 
Z 2 
e it 
ge--¥_(4- + V aed ed )-¥< | (3.5) 
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Seri) 
Ly 1 
(3.6) 
See a er) 
L. 2 
and 
Var ae CE) 
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es ham oe 


Taking the inverse of equations (3.7) the characteristic equations 


of the two nonlinear resistors take the form 


; Sst! 
1 
(3.8) 
; =. yok 3 
2 
where 
1 
h(-) =g. (-) (329) 
Substituting equations (3.1) and (3.2) into (3.5) gives 
ae be els any ue +v ca (3510) 
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Substituting equations (3.1) (3.2), 2@.6)andet3'8) into) (3.3) add 


(3.4) gives 
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Ee Re R 
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Bo 89 
dd dd 
Putting be re and 8 Rivage equations (3.15) and (3.16) take 
the form 
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Examination of the two differential equations (3.17) and (3.18) show 
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that, except for the difference between F,(-) and F,(-), R + R 
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the two equations look the same except that in one half cycle 


R e 
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ot 
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Rt 81 
in the second and this action is reversed in the next half cycle. 
Actually this is the same as saying that in one half cycle the flux 
due to Eo adds to the a.c. flux in the first winding and opposes 
that in the second and the situation is reversed in the next half 
cycle. 

The two differential equations are nonlinear and have 
no analytical solution. In order to use the digital computer to obtain 


a numerical solution it is necessary to formulate these differential 


equations in the normal form [12]. 


a at) (3.19) 


Virtually all numerical techniques for solving a sequence of nonlinear 
differential equations are for those in such normal form. So, in 
order to use any of these numerical techniques it is necessary to 


put equations (3.17) and (3.18) in the normal form. This means 
di dA 


putting an and dp as two explicit functions of hie ro and time 
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da, 
oe E. 4 » Ags t) G2) 
da, dir di, 
ok both F aoe — i i —— 
BS: ro) 9 ie ) and Fy Oe ) are nonlinear functions, dt 


and ar «(Cannot be put in the form of equations (3.20) and (3.21). 
The nonlinearities of FA (+) and F,(+) result from the two nonlinear 
resistors R) and R,- 

As mentioned before in chapter 2 only eight data points 
are available on the characteristics of the nonlinear resistors Ry 
and R,. For calculating intermediate points an interpolation method 
must be used. The piecewise-linearization technique is one of these 
methods. The use of this method, unlike other interpolation 
techniques, makes it possible to put equations (3.15) and (3.16) in 
the required normal form. Applying this piecewise-linearization 
technique the v - i characteristic curve of each of the two non- 
linear resistors is replaced with a number of linear segments with 
the eight data points coinciding with the break points, Fig. 3.5(b). 


Hence, the characteristic equations of the two resistors Ry and R, 


take the form, [APP.A], 
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. = h(v,) = ne? Vy + m2 * pia 12,7 .e, o (323) 
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Substituting these expressions for h(v,) and h(v,) ea Wen CS baal lo Pe 


(3.16) and rearranging give 
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dix ' 
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where 
IT = determinant of the matrix = 


The values of K and m are dependent on the values of nh and ey re- 
spectively. For all possible combinations of m and K, the deter- 
minant IT, has a finite non-zero value. To show that the two equations 
(3.26) and (3.27) have a unique solution that satisfy the initial 


conditions at t = 0, the Cauchy-Lipschitz theorem [14] is used. The 


two functions £(A,) and £0.) are continuous, single valued functions 
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of hy and do respectively. The sinusoidal signal E sinwt is also 
continuous, single valued function for all values of time t. There- 
fore the right hand sides of equations (3.26) and (3.27) are 
continuous, single valued functions of hie do and t for all values 
of Ai hos t and hence satisfy the first condition of the Cauchy- 
Lipschitz theorem. Moreover these two right hand sides are 
differentiable functions of hie do and t. Therefore these two right 
hand side functions are Lipschitz functions and hence satisfy the 
second condition of the Cauchy-Lipschitz theorem. Therefore equations 
(3.26) and (3.27) have a unique solution for all values of t 20. 

D os va and Vv, as functions of time t. 


These two differential equations are solved numerically because there 


This solution provides i 


is no analytical solution available. As mentioned before the object- 
ive is to examine the accuracy of the mathematical model in pre- 

dicting the performance of the parallel-connected magnetic amplifier. 
This can only be done by comparing the actual amplitudes and wave - 
forms of the quantities measured experimentally with the corresponding 
computed amplitudes and waveforms resulting from the mathematical 
analysis in which the mathematical model is used. Therefore, the 
quantities to be compared must be experimentally measurable. But the 
voltages V_ and V, cannot be measured experimentally due to the presence 


dL, 2 
of the finite values of R. and R_. Therefore, the experimentally 
a 3 . 
measurable quantities whose waveforms and amplitudes are to be compared 


with their respective computed waveforms and amplitudes are v4 >the 


voltage across the two parallel-connected gate windings, and the 
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load voltage v,. These voltages are related tov,,v. and e, as 


ee eZ 


BSinwt va V5 

ne = R. SPR sP Pi — Rak (e278) 
1 2) 

Vig Sesingt ety (5229) 


3.3 Numerical computations 


To start numerical computation, using any of the numerical 
techniques to solve the two differential equations (3.26) and (3.27), 
it is necessary to know the initial values of aT and do at t = 0, 
namely, d, (0) and 5 (0)- Before the a.c. voltage source e. is switched 
on at t = 0, there is no a.c. voltage applied to the circuit and 
hence there will be no a.c. current flowing through either one of 
the two gate windings. Under such condition the only flux in the 


two cores will be the d.c. flux produced by the d.c. current flowing 


through the control winding therefore, 
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X, (0) = C, 
The values of the constants Cy and C. are dependent on the currents 
flowing through the two inductors Ly and L, respectively at t=0. 


Since v, (0) = v, (0) = 0 then from the characteristics of the resistors 


R. and R 
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also 


and 


All these imply that the initial currents through the inductors L, 


and L, are 


The characteristics of the two inductors show that the corresponding 


initial flux d, (9) and d, (9) are 


E 
-1 C 
, (0) = dy (0) =f (20 Re 


These initial conditions satisfy the two differential equations at t=0, 
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this can be seen from equations (3.15) and (3.16). 

A computer program for solving equations (3.26) and 
(3.27) was written. This program consists of a main program and 
four subroutines in addition to the library subroutine HPCG [15]. 

The library subroutine HPCG is one of the subroutines 
that can be used to solve a sequence of first order ordinary 
differential equations with given initial values. This subroutine 
uses the Hamming predictor-corrector method. This method is not 
self starting, that is, the functional values at a single previous 
point are not enough to get the functional values ahead. Therefore, 
to obtain the starting values, a special Runge-Kutta procedure 
followed by one iteration step is used added to the predictor-corrector 
method. The Hamming predictor-corrector method is a stable fourth 
- order integration procedure that requires the evaluation of the 
right hand side of the sequence of the differential equations only 
two times per step. This is a great advantage compared with other 
methods of the same order of accuracy, especially the Runge-Kutta 
method, which requires the evaluation of the right hand side four 
times per step. Another advantage is that at each step of calculation 
procedure this subroutine is able, without a significant amount of 
calculation time, to choose and change the step size. 

The subroutine FCT is one of the external subroutines 


that should be used with HPCG and it is intended to compute the right 
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FIG. 34.5 


The i-\ curve (a) for the nonlinear inductors 


and V-i curve (b) for the nonlinear resistors 
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hand side of the two differential equations (3.26) and (3.27) for 
dir di, 
iven value £2 CEMA ie@. sag ete 
g ues oO » Ay and dos ise. calculate ae and FRE 


The subroutine. AVR is used to determine the operating 
point on each of the piecewise-linear characteristics of the two 


resistors R) and R, corresponding to voltages oy and Vv, respectively. 


This is done by determining the segment on which the operating 
point lies for each of the two resistors, thus the values of K, m and 


hence _G 


Cen 


9 corresponding to the operating point 


are known. 
As in all practical cases the characteristics of each 


of Ly and L, are known only by a finite number of experimentally 


obtained data points, in the present work these are eight, and an 
interpolating subroutine should be used. The subroutine BWLINT is 
an interpolating subroutine to enable the calculation of values 


My and do on the characteristics of Lj and L, for given values of i 
1 


and it » and vice versa. This routine consists of only first order 
2 


terms and is equivalent to a piecewise linear representation for the 


characteristics of Ly and Los Fig. 3.5(a), with the data points 


coinciding with the break points [APP.A]. 


The second external subroutine required by HPCG is the 
dr 
1 
output subroutine OUTP. This is intended to store hi No» FES 
dd 


aoe and t as an array and then to plot Mi and ¥e versus time t. 
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The main program and its subroutines are used to solve 
equations (3.26) and (3.27), calculate and plot the steady state 
values of he and v,, versus time. This was done for values of Eo 
equal to 0.0, 5, 10, 15, 20, 25 and 30 volts. The resulting wave- 
forms are shown in Figs. 3.6(a) and (b) to 3.12(a) and (b). The 
a.c. voltage source 2s has an amplitude equal to 17.5 volts and 
frequency of 60 Hz. The values of the resistors used in the 
calculation are exactly equal to those used in the practical circuit 


of the parallel-connected magnetic amplifier. 


3.4 Experimental results 


The values of the parameters of the actual parallel- 


connected magnetic amplifier are 


sme 17.5 sin (27 60t) 
R. = 100 ohm 
Ro = 10189 ohm 
sis 
N. = 20 
& 
R= 0.76 ohm 
ui 
Ro° #05765 ohm 
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The values of the d.c. control voltage is taken as O02 oi eLO..e aie. s 
20, 25 and 30 volts. For each of these values the steady state wave- 
forms for Vi? the voltage across the load R > and Vv the voltage 
across the parallel connection of the two gate windings, were photo- 
graphed. These waveforms are shown in Figs. 3.6(c) and (d) to 


3.12(c) and (d) for each of these the time scale was 2.5 ms /division 


and the voltage scale was 10 volts / division. 
3.5 Comparison 


The checking of the accuracy of the mathematical model 
in predicting the accurate performance of the parallel-connected 
magnetic amplifier is the objective of this section. The waveforms 
calculated and plotted from the analysis, in which the mathematical 
model replaces the magnetic iron-coreg, are compared with the actual 
experimental waveforms.. This was done for values of Eo equal to 0.0, 
5, 10, 15, 20, 25 and 30 volts. The waveforms are shown in Figs. 3.6 
to 3.12 where Figs. (a) and Figs. (b) are the waveforms of the load 


voltagev , and the gate voltage es respectively obtained using the 


L 
mathematical model while Figs. (c) and Figs. (d) are the actual 
experimental load voltage ‘. and gate voltage - waveforms. in-ail 


these figures the voltage scale is 10 volts / division and the time 


scale is 2.5 ms / division. 


Fig. 3.6 was obtained for Eo = 0.0 volts. It is clear 


that the waveforms obtained using the mathematical representation are 
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in good agreement with the experimental waveforms. Actually the 
only difference is the presence of the sharp points appearing in 
Fig. 3.6(a) and (b). The sharp points at the maximum and minimum 
values of the load voltage in Fig. 3.6(a) results in the apparent 
difference in the amplitude of the waveforms in Fig. 3.6(a) and (c). 
Concerning the waveforms it is clear that the mathematical model 
successfully predicts when the gate does or does not take most of 
the a.c. voltage. Fig. 3.7 was obtained for ee = 5 volts. The 
duration over which the load takes most of a.c. voltage starts to 
increase and thus the a.c. power delivered to the load increases. 
Comparison of Fig. 3.7(a) and (b) with Fig. 3.7(c) and (d) re- 
spectively shows that the agreement concerning both amplitudes and 
waveforms is excellent except for the previously stated difference of 
the sharp points appearing in Fig. 3.7(a) and (b). 

As Eo increases, the duration of time over which the load 
takes most of the a.c. voltage start to increase. Figs. 3.8, 3.9, 


3.10, 3.11 and 3.12 are for E, equal to 10, 15, 20, 25 and 30 volts 


C 
respectively. In all these figures it is clear that the mathematical 


model is very successful in predicting the waveforms of vr, and Wes 


Fig. 3.12 was obtained for E, equal to 30 volts. For this value of 


C 
Eos or greater values, the load has all the a.c. voltage across its 
terminal. This is clear in both Fig. 3.12 (a) and (c). The only 


difference is that small sharp dip in the load voltage in Fig. 3.12(a) 
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does not clearly appear in Fig. 3.12(c). The same holds for the 
small sharp pulse in Fig. 3.12(b) which does not clearly appear in 
Fig. 3.12(d). This difference may be due to the piecewise linear- 
ization technique which causes the slight gradual deformation in the 
sinusoidal waveforms in Fig. 3.12(c) to appear as a sudden dip in 
Fig. 3.12(a). Thus, it is clear from the comparison that the 
mathematical model used in the analysis successfully predicts the 
performance of the parallel-connected magnetic amplifier with no 
feedback with the exception of the sharp voltage changes appearing 
in those waveforms computed using the model. This may be arising 
from using the piecewise linearization technique, in interpolating 
? Los Ri and Ry > 


of the model. The exactness of the model may be improved by in- 


the characteristics of L together with some inaccuracy 


serting some parasitic elements. 
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Experimental and computed voltage waveforms for Eo = 0.0 
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Experimental and computed voltage waveforms for Eo = 5 
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Experimental and computed voltage waveforms for Eo 
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Experimental and computed voltage waveforms for Eo = 20 
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Experimental and computed voltage waveforms for Eo = 25 
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CHAPTER 4 


SERIES-CONNECTED MAGNETIC AMPLIFIER 


4.1 Theory of operation 


The circuit diagram of a series-connected non-feedback 
Magnetic amplifier is shown in Fig. 4.1 [10-11-13]. The same two sat- 
urable reactors used for the parallel case is employed. The two 
gate windings are connected in series with their polarities such 
that any a.c. voltage applied across them induces no current in the 
contral circuit. On the load side of the series-connected magnetic 
amplifier the a.c. voltage source e, is in series with the load 
and in series with the two gate windings, the magnetic amplifier 
purour is the a.c. voltage across the load R- If d.c. output 
is required a bridge rectifier can be used. The control side is 
similar to that of the parallel-connected magnetic amplifier with the 
exception that the value of R, is lower. For E, = 0, no d.c. current 


C C 


flows in the control winding N, and each of the gate windings will 


C 
have large inductive impedance. Thus most of the a.c. voltage Ss 
is across the gate windings and a negligible a.c. voltage appears 
across R,- For Eo + O,in one half cycle the flux linkage due to 
the d.c. control current i, adds to the a.c. flux linkage of the 


first gate winding, and opposes that of the second gate winding. 


The high value of flux linkage of the first gate winding drives it 
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towards saturation and the impedance across this first gate becomes 
very small, practically a short circuit. At the same time there is 
a low value of flux in the second gate winding causing it to have 


a high input impedance. Due to the transformer action this high 


N 
a2) 
N 


Z 


impedance is in parallel with the low value of resistance Ra 


resulting in a low impedance across the terminals of the second 

gate winding. Therefore, during this half cycle the impedance 
across the series connection of the two gate windings remains small 
and most of the a.c. voltage appears across the load. During the 
next half cycle the same action takes place but the two gate wind- 
ings replace each other in action and still most of the a.c. voltage 
appears across the load. Thus in the presence of Eo most of the 


a.c. voltage appears across R The amplitude of the load voltage 


together with its waveform are completely dependent on the magnitude 


of Eo: 


4.2 Analysis 

As in the case of the parallel-connected magnetic amplifier, 
the first step in the analysis is to replace each of the two iron- 
core transformers, Fig. 4.1, with the equivalent nonlinear lumped 
circuit model, Fig. 2.4. The characteristics of the nonlinear 
resistors and inductors are the same as for the parallel case since 
these do not depend on the external connection of the transformers. 


Thus, the series-connected magnetic amplifier is completely represented 
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by the equivalent circuit shown in Fig. 4.2. Replacing the ideal 


transformers with the controlled source representation, the 


> 


equivalent circuit takes the form shown in Fig. 4.3. v Vos e 
Ss 


5) 
N 1 


; : ; ; C 
eke Eade ge Bo a i. , and — are still the same as 
ae C Ly L, Ri» R, N 


explained before in the parallel case. 
With reference to Fig. 4.3, applying Kirchhoff current 
law to nodes 1 and 2 the currents in the nonlinear inductances 


are, 


(4.1) 


Hs 
i] 
| 
+ 
No 
oO 
4. 
i=) 
i 


1 
It 
I 

H 

+ 

nN 

=) 

}- 
I 

# 


(4.2) 


With reference to Fig. 4.3, applying Kirchhoff voltage law and 


rearranging the resulting equations, the control current ip and the 


load current i, are 
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(4.4) 


Substitution in equations (4.1) and (4.2) by equations (3. Gi) 


and (4.4) the currents i, and i. are 
1 2 
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(4.5) 


(4.6) 


(4.8) 


(4.9) 
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Equations (4.8) and (4.9) are two nonlinear coupled differential 


equations in ST and Noe As stated before no analytical solution 


is available and in order to solve these two equations numerically 
they should be put in the normal form. As in the parallel case 
this is impossible without using the piecewise linearization 


technique to relate the nonlinear resistors currents in and in 
u Y 
to the voltages %) and V5 respectively. As a result of this 


piecewise linearization technique the two currents ip and ip are 
1 2 


put in the form previously given 
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As before the values of K and m are dependent on the values of v4 and 


V5 respectively. Substitution in equation (4.8) and (4.9) for h(v,); 
h(v,) and collecting terms of v1 together and those of Vo together the 


two equations can be put in the forn, 


E Sint 20E., Vy 
CS CIE "ad Ga a! 
= Ts (4.10) 
E ,sinut 20E,, V5 
£(A,) 2 an ~ap at ane 
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Where 
eons ms 400 in G 1 400 
R Ro Kea: R Ra 
qT, = (ee 
ou, 400 ‘ pron 400 =. 
R Ro R Ro m 2 


Taking the inverse of this transformation matrix, VY, and V, are 


1 2 
given as 
dx E sinwt 
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Where 
IT, | = determinent of the matrix [T_]. 
These two differential equations are in the normal form relating the 


fluxes aT and do with time t. The values of K and m are dependent 


on the values of v1 and v5 respectively. The determinant ee has 


a finite non-zero value for all possible combinations of K and m. 
Similar to the parallel case it is found that the two right hand 
sides of equations (4.12) and (4.13) satisfy the two conditions of 
the Cauchy-Lipschitz theorem, thus these two equations have a unique 


solution for all values of t>0. This solution provides ie os v4 


and V5 as function of time t. But Vi and V5 cannot be measured and 


the measurable quantities in case of the series-connected magnetic amplifier 
are the voltages across the two gate windings we 5 ve > ye and the 


L Z 
load voltage er These are related to the voltages Vi> Vo and e, as 


R R R 
124 ; 81 8} mer 
Vg -E, pg Sin wt ty [1 - rp ie’ aaa (4.14) 
R R R 
ade dy ane = (4.15) 
os. = nie zp sin wt + V5 [1 - rp. ViR : 
v =v -v (4.16) 
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= Eg sinwt - vs 4.17) 


4.3 Numerical Computation 


To start the numerical computation, it is necessary to calculate 


the initial values of aT and ho at t = 0, namely, , (0) and 5 (Q)- 


Following the same reasoning as in the case of parallel-connected 


Magnetic amplifier it is found that, 


v, (0) = vy (0) = 0 (4.18) 


OG, a ee (26S) 
160 = (0 Ro (4.19) 


These initial conditions clearly satisfy the equation (4.10). 
The same computer program used in the case of the parallel-connected 
Magnetic amplifier is used here with the exception that the subroutines 


FCT and OUTP are changed. The new FCT subroutine is written to 
da dd, 
compute the values of 7 and ae for any given Aas h, and t using 


eguations (4.12) and (4.13). The new OUTP subroutine calculates and 


stores values of the load voltage v,, the gate voltages VE Gave ana Mee 
= it & D 


then plots v hg and y_ versus time t. 


_ 
The steady state waveforms of v,, v, and She were plotted versus 
i a 
time for values of E, equal to 0.0, 10, 20, 30 and 40 volts. These 


[ai 


waveforms are shown in Figs. 4.4(a), (b) and (c) to Figs. 4.8(a), (b) 
and (oe). oneea.c.,. voltage source C4 has an amplitude equal to 25 volts 


and frequency of 60 Hz. The values of the resistors used in the 
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calculation are exactly equal to those used in the practical circuit 
of the series-connected magnetic amplifier. 


4.4 Experimental results 


The values of the parameters of the circuit used in experimental 
investigation are the same as those used in the parallel-connected 
Magnetic amplifier experimental investigation except that the value 
of Ra now is 5 Kf and amplitude of e is 25 volts. 
The values of the d.c. control voltage are taken as 0.0, 10, 
20, 30 and 40 volts. For each of these values the steady state waveforms of 
vy? Ne and v : were photographed. These waveforms are shown in 
Rigel: 4.4(d), (e) and (f) to Figs. 4.8(d), (e) and (f). In each 
of these photos the time scale was 2.5 ms /division and the voltage 


scale was 10 volts/division. 


4.5 Comparison 


The accuracy of the mathematical model in predicting the accurate 
performance of the series-connected magnetic amplifier is checked by 
comparing the waveforms computed and plotted from the analysis with 
the actual experimental waveforms. The waveforms of three voltages 
are compared, namely, the load voltage Vi? the voltage across the 
terminals of the first gate winding Va and the voltage across the 


two gate windings v_. The comparison was carried out for values 
& 


of the dsc. control voltage Eo equal to 0.0, 10, 20, 30 and 40 volts. 
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The waveforms corresponding to these values are shwon in Figs. 4.4 to 4.8. 
Figures (a), (b) and (e) show the waveforms of Vi o and vs re- 
spectively, computed using the mathematical model. Figures Gn (e) 

and (f) show the actual experimental waveforms of v., ne and ae wee ot 

is clear that the amplitude of tay is nearly half that of ve a 

is also clear that for either the experimental waveforms or those 


is always 


plotted from the modeling analysis, the sum of rs and vs, 


equal to the sinusoidal voltage Ey sin (27 60)t. 

Fig. 4.4 was obtained for E,=0.0 volts. The load voltage is 
of a very small amplitude while the gate voltages are nearly sinu- 
soidal. The agreement concerning both amplitudes and waveforms is 
excellent except for the sudden sharp changes in the waveforms in 
Pigs. -4.4{(a), (b):-and (c). 

As the magnitude of the d.c. control voltage Eo increases, the 
load resistor sae starts to have most of the a.c. voltage across its 
terminals and the duration over which this occurs increases as Ee 
increases, Figs. 4.4 to 4.8, till the load voltage v, becomes 
nearly sinusoidal at Eo = 40 volts, Fig. 4.8. Figures 4.4 to 4.8 
show excellent agreement between the waveforms of v_, ve and v ‘ in 
Figs. (a), (b) and (c) and those in Figs. (d), (e) and (f) respectively. 
The agreement is excellent concerning the amplitudes of these voltages 


as well as their waveforms. 


Thus except for this slight difference, in the form of sudden 
sharp change of the voltage in those waveforms computed using the model, 


the mathematical model successfully predicts the accurate performance. 


8.8 ot &.\ Vagkt ot nowde ox’ 


<97 ,.¥ bab V @4¥ to me ots Ay Bie ad wore (3). brs ‘the 


Pas | | aes 
(os) .€by e5tyait. {ebom {dot semod sea aid. gareu besuqnos, e 


40) OV) ies? UR oe 50. aim Yevew JIeinsmbreqxs, teutos odd. word» ; 


x ; S A i" 


4¥ . ¥ to tadd tish wisasit at igvdo akatis queg Poe eda 
> q ‘a 
} 2 J ; ae f cet f 


sa0ds xtolamressvew [s3comh7 gus adt sedate ee sie. xaos 

5 | as = te : : Ay) a Ly 

gyswis ak .v bos /V/40) mye ods ereuisas! gaklebom eri cont. 
‘ | ie ‘ey * aki 


ax fk . ae yp 
<3(0d nS), nate 3 oasiidy ‘Tabioeunbe add oF 
at Bee rLov beol si?’ .etted Ore ¢ 402 dhonbes ae asw' As) ght 

; ee MS, ge Riza oe 
veunte vitser s1B eseygstioy sites efi siitiw sbutiigms fisme yxsv) 


ei eartotevew bre sebus i lqas aod galerie smo; Hema SE: 


nt amrotsyvaw sii a asymp qrsde mabbua adi x02 aqaowe fe 7“ 


=u 


lag é. ee edad 4 


a “odd -esasoront .o) egsd.loy Loxtnoo. ‘Died oda qo sbusingam ads 
| ae Tea hy i | 


 -gdk eeotos egstiov ios. aa to. Jeon ove! or er see rodetess bao! 
| | | “1 * a 


tT ae sésazant a0990 aida do kite. asvo’ notgexub eda ein ‘elsah 


vasmosed yr) sgsatov | beot ous reer 3. i od Aa s .aget — 


| J a, oe = 
Bh of vA eozuatt 8. & att eitoy on er ae Aebtosuate | in 


; gk gr aM et vA he To aarrolavew oda moeyiad snomoorae sats: 

ri : 7” eo ae " "i 

pvisoedeos Q) brs on by att ak spot bas. ee bas @o 
oo) Ln 


- i Dy gi an fo he + 
; HL: Teer bee st bie ee es 


[J 


7 4 ; faba ye gates bonugao amotavus ovo Ps 8 reat is . 
Carls * a ei 

J A as . A a 7 7 q Pu 
Pts ae sseti208 eas ssolbosa vehotaneee ue Lebom Ls 


i a n ' 7 Sa wl " . a os a 


Ah aoe RD, 


67 


of the series-connected non-feedback magnetic amplifier. As stated 
in the parallel case this may be the results of using the piecewise 
linearization technique in interpolating the characteristics of the 


nonlinear elements together with some inaccuracy of the model. 
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FIG. 4.4 
Experimental and computed voltage waveforms for E, = 0.0 
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FIG. 4.6 


Experimental and computed voltage waveforms for Eo = 20 
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Experimental and computed voltage waveforms for Eo = 30 
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Experimental and computed voltage waveforms for Eo = 40 
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CHAPTER 5 
CONCLUSION 

The mathematical model used in this work to represent the 
magnetic iron-cores was proved once to successfully predict the 
correct performance of a push-pull tunnel-diode relaxation oscillator 
with a single iron-core. This model was found to be a very simple 
and useful one especially in representing magnetic iron-cores since 
it yields a lumped circuit model. Since this is a postulated non- 
linear mathematical model, it is necessary to prove that it can 
predict the accurate performance of a large number of systems comprising 
hysteretic elements. In the present work it is found that the model 
successfully predicts the accurate performance of both the parallel 
and the series connected magnetic amplifiers. 

First the parameters of the lumped circuit model were obtained 
from the hysteresis loop of the magnetic iron-core. Then the lumped 
circuit model replaces the two magnetic iron-cores. The resulting 
circuits for the parallel and the series connected magnetic amplifiers 
were analyzed and in both cases the resulting two differential 
equations were solved numerically. 

The actual experimental waveforms were then compared to those 
computed from the analysis and it was found that for both types of 


connections the model is very successfull in predicting the experimental 
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waveforms except for a very slight difference in the form of sharp 
points in the computed waveforms, which may be due to the use of 

the piecewise-linearization techniques in interpolating the 
characteristics of the nonlinear elements of the model together with 
any slight inexactness that the model may have. This slight in- 


accuracy may be avoided by adding some parasitic element to the model. 


The only drawback in the simulation process of the parallel- 

connected and the series-connected magnetic amplifier circuits is the 

cost of the numerical computations. In case of the parallel-connected 
magnetic amplifier, steady state is reached after about 45 minutes of 

epi cation on the IBM 360 digital computer. In case of the series- 
connected magnetic amplifier the required time, of course smaller than 

in the parallel-connected case,is about 5 minutes. The time of computation 
may be reduced by finding more suitable initial conditions and using 


a numerical technique faster than the HPCG. 
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APPENDIX A 
THE PIECEWISE LINEARIZATION 
TECHNIQUE 
If only a finite number of points on the characteristics of 
a nonlinear element are known and it is required to calculate the 
coordinates of intermediate points, it is necessary to use some 
interpolating technique [12]. The piecewise linearization tech- 
nique replaces the actual nonlinear curve, which join this point, 
with a number of linear segments. The break points of these segments 
are the known data points. 
If the two variables are x and y and the number of given 
data points is ee then the nonlinear curve is replaced by So. 
segments, where j is the nonlinear element number. Each of the 
segments is completely specified by: 
1\— Los of the asa segment = Ke; 
2 - Intersection of the nes segment with the y axis = Ky 
3 - Interval of definition of the ea segment, namely, 
rm 2 + 


Yenc acy eee Ne. Yonik ee 
Ko RR he? WK ae Bs 


Fig. (A.1) shows a typical ron segment illustrating various symbols 
defined above. 

If the nonlinear element is a resistor the variables x and y are 
the voltage VR and the current ips If the nonlinear element is an 


inductor the x and y are the current i and flux i. 
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The parameters defining a piecewise-linear segment 
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APPENDIX B 


FLOW CHARTS 


MAIN PROGRAM 


START 


READ NUMBER OF BREAK 
POINTS FOR NONLINEAR 
RESISTORS AND INDUCTORS 


READ AND WRITE COORDINATES OF 
BREAK POINTS FOR NONLINEAR 
RESISTORS AND INDUCTORS 


READ THE INITIAL INCREMENT OF 
TIME (DELT), UPPER ERROR BOUND 
(ACC), NUMBER OF Ec VALUES 
(MAXE) AND MINIMUM TIME (TM) 


READ AND WRITE 
MAGNETIC AMPLIFIER 
CIRCUIT PARAMETER 


CALCULATE D.C. CONTROL CURRENT 


CALCULATE AND WRITE THE CONDUCTANCE AND 
INTERCEPT WITH THE CURRENT AXIS FOR EACH 
SEGMENT OF THE PIECEWISE LINEARIZED 
RESISTORS CHARACTERISTICS 
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PUT INITIAL VALUES AT T=0.0 
T=0.0, V(1)=0.0, V(2)=0.0AND 
i, (1) =i, (2) =CONTROL CURRENT 


CALL BWLINT 
"CALCULATE INITIAL 
INDUCTORS FLUX" 


WRITE INITIAL 
INDUCTORS CURRENTS 
AND FLUXES 


SET PARAMETERS REQUIRED BY 
THE LIBRARY SUBROUTINE HPCG 
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SUBROUTINE BWLINT 


START 


CALCULATE FLUXES CORRESPONDING TO GIVEN 
) VALUE OF INDUCTORS CURRENT AND VICE VERSA 


GIVEN 
CURRENT (FLUX) OF INDUCTOR (II) 
SMALLEST -VE CURRENT (FLUX) OF BREA 
POINTS OF INDUCTOR (II 


GIVEN 
CURRENT (FLUX) OF INDUCTOR (IT) 
LARGEST +VE CURRENT (FLUX) OF BREAK 
POINTS OF INDUCTOR (II 


CALCULATE FLUX (CURRENT) 
USING THE EQUATION OF 
SEGMENT (1-1) 


FLUX (CURRENT) = FLUX (CURRENT) 
AT BREAK POINT (1) 


YES 


RETURN 
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SUBROUTINE FCT 


CALCULATE THE ELEMENTS 
OF THE IMPEDANCE MATRIX 


CALL BWLINT 
"CALCULATE INDUCTORS 
CURRENT“ 


CALCULATE THE ELEMENTS 
OF THE CURRENT COLUMN 
CALCULATE VOLTAGES 
V(1) AND V(2) 


RETURN 
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SUBROUTINE AVR 


START 


CALCULATE THE PARAMETERS FOR THE 
TWO SEGMENTS, OF THE RESISTORS, 
j USED ATASPECIFIC INSTANT 


VOLTAGE OF 
RESISTOR (IA) < THE SMALLEST 
-VE VOLTAGE OF BREAK POINTS 
OF RESISTOR (IA 


VOLTAGE OF — 
RESISTOR (IA) > THE LARGEST 
+VE VOLTAGE OF BREAK POINTS 
OF RESISTOR (IA 


S VOLTAGE 
OF BREAK POINT (MM) MINUS 
VOLTAGE OF RESISTOR (IA 
-,0, as 


CALCULATE SLOPE AND 
CURRENT AXIS INTERCEPTION 
OF SEGMENT (MM-1) 


RETURN 


1) «edie l can Sr 


oa 4 


SUBROUTINE OUTP 


TMAX = TM+.025 
NO at 
¥ 


ES 
MAAN = MAAN+1 


CALCULATE INSTANTANEOUS 
VALUE OF A.C. VOLTAGE 


STORE TIME, FLUXES, 
VOLTAGES V(1) AND V(2) 


CALCULATE AND STORE 
GATE AND LOAD VOLTAGES 


MAAN1 = MAAN+1 
MAAN2 = MAAN+2 


“$_<> 
YES 


PLOT THE REQUIRED 
VARIOUS GRAPHS USING 
CALCOMP PLOT SUBROUTINES 


RETURN 
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